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Abstract 

In the semiclassical approach, inclusive and diffractive quark and gluon distribu¬ 
tions are expressed in terms of correlation functions of Wilson loops. Each Wilson 
loop integrates the colour field strength in the area between the trajectories of 
two fast partons penetrating the proton. We introduce a specific model for aver¬ 
aging over the relevant colour field configurations. Within this model, all parton 
distributions at some low scale Qq are given in terms of three parameters. Inclu¬ 
sive and diffractive structure functions at higher values of are determined in a 
leading-order QCD analysis. In both cases, the evolution is driven by a large gluon 
distribution. A satisfactory description of the structure functions F2{x,Q‘^) and 
S,Q‘^) is obtained. The observed rise of with ^ is parametrized by 

a non-perturbative logarithmic energy dependence, compatible with unitarity. In 
our analysis, the observed rise of F 2 at small x is largely due to the same effect. 


1 Introduction 


Recently, precise measurements of inclusive |l| and diffractive p| structure functions 
at small x have become available. Numerical analyses, based on different theoretical 
approaches, have been performed by many authors (see e.g. p-§[). 

In the present paper, a joint analysis of diffractive and inclusive structure functions 
based on the semiclassical calculations of [|^,H] is performed. The idea of a close similarity 
between diffractive and non-diffractive processes in deep inelastic scattering (DIS) lies 
at the heart of the semiclassical approach. In both cases, a partonic fluctuation of the 
incoming virtual photon scatters off a superposition of target colour helds. If the scattered 
partons emerge in an overall colour singlet conhguration, a diffractive hnal state results. 


Comparing the semiclassical description of structure functions with parton model 
expressions, we dehne inclusive and diffractive parton distributions in the semiclassical 
approach. Higher order contributions in the semiclassical calculation exactly reproduce 
the leading logarithmic corrections to the parton model, showing the consistency of both 
approaches. The semiclassical method is therefore used at some low scale to derive initial 
distributions. Starting from these distributions, a leading-order DGLAP analysis [0 of 
experimental data is performed. 


The calculation of the above initial distributions involves averaging over all relevant 
colour held conhgurations. To perform this averaging, a simple non-perturbative model, 
valid in the case of large hadronic targets [M, is used. It is based on the observation 


that, for extended target colour helds, the transverse size of partonic huctuations of the 
photon remains small . 


The energy dependence arising from the large-momentum cutoh applied in the pro¬ 
cess of colour held averaging can not be calculated from hrst principles. It is described 
by a In a: ansatz, consistent with unitarity. In the semiclassical approach, this energy 
dependence is expected to be universal for both the inclusive and dihractive structure 
functions |T^. 


Overall, a satisfactory description of inclusive and dihractive small-x structure func¬ 
tions, based on a minimal number of assumptions and only four htted parameters, is 
achieved. In our opinion, this lends support to the idea of a close similarity between 
the mechanisms of inclusive and dihractive scattering and its implementation in the 
semiclassical approach. 


The paper is organized as follows. In Sections ^ and |^, semiclassical formulae for 
inclusive and dihractive parton distributions are given, and the underlying physical pic¬ 
ture is discussed. Section ^ deals with our model for the colour held averaging that is 
responsible for the input distributions. The leading-order DGLAP analysis and the com¬ 
parison with experimental data are the subject of Sect. |^, followed by conclusions in 
Sect. p. Appendix A contains some additional formulae relevant for the calculation of 
the dihractive parton distribution functions. Finally, we illustrate in Appendix B how 
the dihractive parton distributions of a small colour dipole, calculated in 0, can be 
obtained using the semiclassical method. 
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Figure |^: Inclusive DIS in the proton rest frame (left) and the Breit frame (right); 
asymmetric fluctuations correspond to quark scattering (a), symmetric fluctuations to 
boson-gluon fusion (b). 


2 Inclusive structure functions 


Small-x DIS can be conveniently discussed in terms of the qq wave function of the virtual 
photon (see e.g. [Q]). In the semiclassical approach, the corresponding qq states scatter 
off a ‘soft’ target colour held (cf. the l.h. side of Fig. |I|). As a; —>• 0, the resulting inclusive 
structure function Ft{x,Q'^) approaches a constant |^, 

Ft(x,Q^) = da(a^ + {l-o:f)N^ J^^K,{yNf 

+ 0(x). (1) 

The light-like paths of quark and antiquark penetrate the colour held of the proton at 
transverse positions x± and x± + y± picking up non-Abelian phase factors U^{x±) and 
U^^{xj_ + y±) (where T stands for the fundamental representation). The function 

= U^{x^)\J^\x^ + l/±) - 1 (2) 

is essentially a closed Wilson loop through the corresponding section of the proton, which 
measures an integral of the proton colour held strength. Furthermore, a is the fraction 
of the photon momentum carried by the quark, = a{l — a)Q‘^, y = \y±\ and is the 
quark charge. 

To map this calculation onto the conventional parton model framework, extract the 
leading-twist contribution from Eq. m and identify it with Ft{x,Q‘^) = 2egXq{x,Q^). 
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The resulting quark distribution reads 


xq{x, Q^) 


+5(|;^(ln^-l)/^tr(5,,lV;(0H,W^;(0)) . (3) 


Note that this expression is independent of which has only been introduced to separate 
the logarithmic dependence explicitly. We have assumed that the hadronic structure, 
encoded in the y'^ dependence of , is dominated by the soft scale A^, which satishes 
< /i^ < QA 

The corresponding gluon distribution at small x is most easily calculated as 


xg{x,Q‘^) 


Stt dFxix, Q'^) 

agCg dlnQ"^ 




(4) 


Equations @ and (|) provide the basis for our analysis of inclusive DIS. Note that 
the semiclassical approach predicts gi^x^Q"^) ~ l/o:, which is expected for a classical 
bremsstrahlung spectrum of gluons. 

To gain more physical insight into the correspondence of the semiclassical and the 
parton model approach, return to our starting point, Eq. (0). It is instructive to view 
Ft as a sum of two terms: the contribution of asymmetric configurations where 

quark or antiquark are slow, a < or 1 — a < (Fig- |l|a), and the 

contribution of symmetric configurations where both quark and antiquark are fast, a, 1 — 
a > (Fig. |l|b). In an infinite momentum frame of the proton, the asymmetric and 

symmetric contribution to Ft correspond to photon-quark scattering and photon-gluon 
fusion respectively. 

The symmetric part is dominated by small qq pairs, i.e., by the short distance con¬ 
tribution to the Wilson-loop trace. 



= tr (9,,Wp0)B,^W^l(0)) + 0(v‘) . 


(5) 


The corresponding contribution to the structure function is related to the second term on 
the r.h. side of Eq. (Q), which generates the gluon distribution, Eq. (|). It was evaluated 
in 0] at leading order in A^//i^ and 

Fr‘[a,Q^) = -§-,jp^P„(z)i\xi^-'^ j^^t,:{d„^wpo)d,^wpya)). (e) 


Here Pqg{z) is the conventional gluon-quark splitting function. 

The other splitting functions appear if as corrections to Ft and Ft, associated with 
higher Fock states of the virtual photon, are considered in the semiclassical approach. 
For example, the qqg parton conhguration involves, in the case where one of the quarks 
carries a small fraction of the photon momentum, a InQ^ term associated with Pqq{z). 
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The splitting function Pgg{z) is most easily derived by considering an incoming virtual 
scalar which couples directly to the gluonic action term . Its lowest order Fock 

state consists of two gluons. We have checked explicitly that the semiclassical calculation 
of the corresponding high energy scattering process yields the usual gluon-gluon splitting 
function. 


Having shown that the leading logarithmic QCD corrections to the semiclassical ap¬ 
proach exactly reproduce the well-known DGLAP splitting functions |]^, we do not 
pursue the calculation of higher order as contributions along the lines discussed above. 
Instead, the large logarithms In (Q^//i^) are resummed in the conventional way, by means 
of the renormalization group. To this end, the parton distributions q{x, and g{x, 
are evaluated using DGLAP evolution equations, with the input distributions q{x, Ql) 
and g{x,Ql) given by Eqs. (|]) and (^). Here Ql is some small scale where logarith¬ 
mic corrections are not yet important. The parton model description of the structure 
function at leading order includes only photon-quark scattering. The leading logarithmic 
term from the photon-gluon fusion process appears now as part of the resummed quark 
distribution. 


So far, we have considered electroproduction off a hxed ‘soft’ colour held. As a con¬ 
sequence, the unevolved structure function Ft{x, Ql) approaches a constant value as 
a; —>■ 0. However, a proper treatment of the target requires the integration over all rele¬ 
vant colour held conhgurations. 

The qualitative features of the held averaging procedure are most conveniently dis¬ 
cussed for the simple case of qq pair production in the proton rest frame. Gonsider the 
corresponding amplitude in a ‘mixed’ representation, < qqA\'j*p >, where the hnal 
state consists of the outgoing qq pair and a colour held conhguration A. We neglect any 
time evolution of the held between the actual scattering process and the moment at 
which the hnal state held conhguration A is dehned. The squared amplitude, summed 
over all helds A and normalized to the total space-time volume, reads 


1 

vf 



< qq A \ 'j*p > 


4:7imp6{k^ + 



%[Aioc\ F[Aioc] 


(7) 


Here <I>p is the proton wave functional, F is dehned by the amplitude for the scattering 
oh a hxed held A, 

< >A= 27r(5(A:J + /c? - g°) F[A\ , (8) 

rrip is the proton mass, and g, kq, kq are the momenta of the incoming photon and the 
outgoing quark and antiquark respectively. The index ‘Zoc’ symbolizes that, on the r.h. 
side of Eq. (|^), the integration is restricted to helds localized at, say, T = 0. This can 
be justihed using translation covariance of the proton wave functional and of the matrix 
element < qq\Y (cf- Sect. 2 of [0). 

When writing IF^(|/j_) we have, until now, always assumed the functional depen¬ 
dence on the classical colour held conhguration A^ to be implicit, so that one should 
really read HA^(|/_l)[Acz]. As can be seen from Eq. (|^), the full inclusive parton distribu- 
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tions are obtained from the previous formulae by the substitution 


trfW'fysJIA,] iuy(!/l)[^d]) - / I MAioJ f tr(wfy9x)[dl,J . 

^ -^loc 

(9) 

The same applies to the diffractive distributions of the next section. An explicit model 
for the held average will be described in Sect. I]. 

Decomposing the held Aioc in Eq. (|^ into its Fourier modes Aioc{k), the path integral 
can be written as 

/ = n jdAUk), (10) 


where the cutoh \q\ is required to ensure that the basic precondition for the semiclassical 
treatment, the softness of the target colour held with respect to the momenta of the 
fast particles, is respected. This cutoh induces a non-trivial energy dependence of the 
squared amplitude in Eq. ( 0 )- 

We are not able to derive the explicit form of that energy dependence from hrst 
principles. Nevertheless, based on the above qualitative discussion, we ascribe a soft, 
non-perturbative energy growth to our input parton distributions used in the numerical 
analysis of Sect. 0. 


3 Diffractive structure functions 


The dihractive cross sections obtained in the semiclassical approach can be expressed 
as convolution of the ordinary partonic cross sections and dihractive parton distributions, 
calculated in [Q. The qq and qqg conhgurations then yield 


da' 


T,L 


di 


j dy\^ + dT,L{y) 

™ dg{x/y,0 
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dq{x/y,i) 

di 


+ drAv) 
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( 11 ) 


Here ^ is related to the dihractive mass M by ^ = xp = x/(3 and f3 = /{Q"^ + M^). 

The corresponding structure functions read 




dq{h,i) 

di 


a 




+ — Pqg{.z) In — + . . . 


27r 




dg{h,i) 

di 


( 12 ) 


a o 2 f^dbj dq{h,i) ^dg{b,^)\ 

f.K.ftQ) = + (13) 


where z = j3/b, and Cp and Tp are the usual colour factors. The physical interpretation 
of the dihractive parton distributions in the Breit frame is analogous to the interpretation 
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of the inclusive distributions. The function df{b,^)/d^ (where f = q, g) is a. conditional 
probability distribution. It describes the probability of hnding a parton /, carrying a 
fraction of the proton momentum, inside a colour-neutral cluster of partons that 
carries in total a fraction ^ of the proton momentum. 


From Eqs. (0) and (|^, it is obvious that the diffractive structure functions satisfy 
the ordinary DGLAP evolution equations (cf. the more general discussion of [^]). The 
perturbative evolution takes place in the variables (3 and Q^] ^ acts merely as a parameter. 
The physical reason for this is intuitively clear: for an arbitrary DIS event the invariant 
hadronic mass is W, and the quark which couples to the virtual photon can be radiated 
by a parton whose fraction of the proton momentum varies from 1 to a; = Q '^/-|- PF^). 
In a diffractive event, a colour-neutral cluster of ‘wee partons’ is stripped off the proton. 
The invariant mass of this cluster and the virtual photon is M. Hence, W is replaced 
by M, and the quark which couples to the photon can be radiated by a parton whose 
fraction of the cluster momentum varies from 1 to /3 = Q'^/{Q‘^ + M^). 


The separation of the colour-neutral cluster of ‘wee partons’ is non-perturbative and 
independent of the perturbative evolution. It is, however, incorrect to visualize a two- 
step process where the colour-neutral cluster is first emitted by the proton and then 
probed by the virtual photon. If this was the case, the two-gluon or two-quark cluster 
relevant in our calculation (cf. the r.h. side of Fig. would necessarily lead to parton 
distributions symmetric in (3 and 1 — j3. A counter example to this is provided by the 
model distributions derived in Sect. 

The diffractive quark and gluon distributions have been determined in |^. In terms 
of Wilson loops in coordinate space, the quark distribution can be expressed as follows 
(cf. Appendix A), 


^ 2b r r y±y± 

d^ ea-brJ {2'KfNjy,,y'^ yy' 


xKdyN)KdyN)[ trWfj2/x)trlFf;(2/;), (14) 

Jx± 

where N^. is the number of colours and 

. (15) 

1 — 0 


The corresponding expression for the diffractive gluon distribution is very similar. Be¬ 
cause of the different colour and spin of the gluon, the Wilson loop is now in the adjoint 
representation, the tensor structure is different and the modihed Bessel function K 2 
appears^. 


dg{b,d) 

di 


d^k\ k'^ 


e(l - bf J (27r)6iV2 




xK2{yN)K2{y'N) [ , (16) 

Jx^ 

^To simplify the colour algebra of Sect. ^ we use the large W limit throughout this paper, with the 
exception of the standard DGLAP evolution in Sect. ^ 
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Figure Diffractive DIS in the proton rest frame (left) and the Breit frame (right); 
asymmetric quark fluctuations correspond to diffractive quark scattering, asymmetric 
gluon fluctuations to diffractive boson-gluon fusion. 


where 

tp(2/±)=^*,-2^. (17) 

yl 

The physical content of both expressions is rather transparent. The modihed Bessel 
functions contain the kinematic effects due to the propagators of the partons penetrating 
the colour field and the Wilson loops represent the interaction with the proton. Figure ^ 
illustrates the correspondence between the semiclassical and parton model view of the 
leading-order processes testing diffractive quark and gluon distributions. 

To illustrate the close similarity of diffractive and inclusive scattering in the semiclas¬ 
sical approach, it is instructive to evaluate the diffractive contribution to the inclusive 
quark distribution. To achieve this, one sets x = and integrates over ^ keeping x hxed, 

QdKx) = - — -. (18) 

After neglecting terms 0(x) and exchanging the integration variable ^ for iV^, the k'j_ 
integration can be performed trivially. As a consequence, the integrand is evaluated at 
y_L = y'j_, giving the simple result 

/ AfW / KfyNf ( trwyfexjtrwytfei). (19) 

[Z7l) J Jy^ Jx±^ 

This expression can also be obtained from Eq. (H) by the substitution 

tr (irJyyJHyyj/x)) - . (20) 
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After this substitution, the integration in Eq. (^) becomes UV-£nite, so that the upper 
limit fj? can be dropped; the InQ^ term disappears since ti = 0. The above 
relation between the inclusive and diffractive quark distribution illustrates most clearly 
the basic idea of the semiclassical approach: diffractive and non-diffractive events are 
both induced by the interaction of the soft proton colour held with a partonic huctuation 
of the incoming virtual photon. The different hnal states are realized by projecting the 
outgoing partonic system onto different colour conhgurations. 


Let us finally note that the semiclassical formulae for diffractive parton distributions 
Pi used in this paper are sufficiently general to accommodate different models of the 
hadron colour held. For example, they can serve as an alternative starting point for the 
derivation of the perturbatively generated distributions of [113, where a heavy quark- 


antiquark huctuation of a photon was used as a model for the target. Taking the colour 
held responsible for trWtrW^ in Eqs. (|I^) and (|I^) to be the perturbative held of a small 
dipole conhguration, the results of are exactly reproduced (cf. Appendix B). However, 
the analysis of the present paper is based on the fundamentally non-perturbative model 
described in the next section. 


4 A Model for the gluon field averaging 


The averaging over the gluon held conhgurations of the target discussed in the two 
previous sections is a complicated operation depending on the full details of the non- 
perturbative hadronic state. However, in the special case of a very large target, a quan¬ 
titative treatment becomes possible under minimal additional assumptions. 


McLerran and Venugopalan have observed that the large size of a hadronic target, 
realized, e.g., in an extremely heavy nucleus, introduces a new hard scale into the process 
of DIS 0]. From the target rest frame point of view, this means that the typical trans¬ 
verse size of the partonic huctuations of the virtual photon remains perturbative [Q, 
thus justifying the omission of higher Fock states in the semiclassical calculation. Note 
that this does not imply a complete reduction to perturbation theory since the long 
distance which the partonic huctuation travels in the target compensates for its small 
transverse size, thus requiring the eikonalization of gluon exchange. 


Within this framework, it is natural to introduce the additional assumption that the 
gluonic helds encountered by the partonic probe in distant regions of the target are not 

[0 


correlated (cf. and the somewhat simplihed discussion in 


Thus, one arrives at 


the situation depicted in Fig. ^ where a colour dipole passes a large number of regions, 
each one of size ~ 1/A, with mutually uncorrelated colour helds Ai ... An. 


Consider the fundamental quantity Wx^{yA.)ij[A] which, after specify¬ 

ing the required representation and appropriately contracting the colour indices ijkl, en¬ 
ters the formulae for inclusive and dihractive parton distributions. According to Eq. (H), 
this quantity is the sum of four terms, the most complicated of which involves four U 
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Figure Colour dipole travelling through a large hadronic target. 


matrices, 

{u^AA] uA^aa]}^^ ( 21 ) 

= {u.AAA---u.AAA ul^^,AAA---uL+,AAA}^^ 

X {C1..+.1 KJ• • ■ Ixiil uA Ixiil ---uA [xi„l},,. 

The crucial assumption that the helds in regions 1 ... u are uncorrelated is implemented 
by writing the integral over all held conhgurations as 



i.e., as a product of independent integrals. Here the appropriate weighting provided by 
the target wave functional is implicit in the symbol /^. 

Under the integration specihed by Eq. (^), the U matrices on the r.h. side of Eq. m 
can be rearranged to give the result 

IJu^AA] ^^±+y± ( 23 ) 

X C/lJAJ ■ • A,^+yJAl H|21il},, ■ 

To see this, observe that the Ai integration acts on the integrand {Uxj^ [Hi] [Hi]}jq/ 

[Hi] [Hijjfc/;/ transforming it into an invariant colour tensor with the in¬ 
dices i'j'k'V. The neighbouring matrices Uxx_ [H 2 ] and [H 2 ] can now be commuted 
through this tensor structure in such a way that the expression {t 4 j^[H 2 ] 

[H 2 ] UA[A‘^}k"i" emerges. Subsequently, the H 2 integration transforms this ex¬ 
pression into an invariant tensor with indices Repeating this argument, one 

eventually arrives at the structure displayed on the r.h. side of Eq. (|^) . 

To evaluate Eq. (p3|) further, observe that it represents a contraction of n identical 
tensors 

Fijkl= f {UxxArn\Ulj_+yA-^ni\}ij {Uxj_+yA^AUAAA}kl , (24) 

A A.m 

where the index m refers to any one of the regions 1 ... n into which the target is 
subdivided. At this point, we make use of the fact that for a sufficiently large target the 
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transverse separations y±_ and y'^ are always small [0. In fact, for a target of geometrical 
size n/A (where n 3> 1 ), the relevant transverse distances are bonnded by y^ 


y '2 _ 


1/nA^. 


Assnming that size and x±_ dependence of typical field confignrations Am are charac¬ 
terized by the scale A, it follows that the prodncts U^JAl^+y^ and Ux^+y'J^l^ close 
to nnit matrices for all relevant y± and y'^. Therefore, it is justified to write 


Ul^+yAM = exp{zT“r(a;±,|/±)[A^]}, (25) 


where T“ are the conventional group generators and /“ are functions of x± and y± and 
functionals of Am- Equation ( |25|) and its y'^ analogue are expanded around yx = y'± = ^ 
(which corresponds to f'^{x±,0) = 0) ^md inserted into Eq. (|^) . At leading non-trivial 
order, the result reads 

F.jti = S.,hi (l - \lCR(y^ + !,“)) + 'r(vW±)TtjTM . (26) 


where Cr is the Casimir number of the relevant representation [Cr = Cf,a) and 7 is 
dehned by 

/ yA)f\xF. y'A = + 0{yy^). (27) 

A A 

Note that the absence of terms linear in /“ and the simple structure on the r.h. side of 
Eq. (^^ are enforced by colour covariance and transverse space covariance. The absence 
of an explicit x_l dependence is a consequence of the homogeneity that we assume to hold 
over the large transverse size of the target. Neglecting boundary effects, we can account 
for the xf integration by multiplying the final result with a parameter hi ~ T? j that 
characterizes the geometrical cross section of the target. 

Substituting the n tensors Fij^i on the r.h. side of Eq. (|^) by the expression given in 
Eq. ( pSP and contracting the colour indices as appropriate for the inclusive and diffractive 
case respectively, one obtains, in the large-A^c limit. 


where (Ir is the dimension of the representation. 


1 

1 


\lCR(yL - v'±? 

+ y'l) 


(28) 

(29) 


Since n is assumed to be large and the typical values of y"^ and y'"^ do not exceed 
1/nA^, the formula (1 —x/n)" ~ exp[—x] can be applied to the r.h. sides of Eqs. (p 8 |) and 
(P^ . Furthermore, contributions proportional to {U^JjA+yAiAku ^ij{Ux^+y'JAl.A}ki 
and 5ij5ki have to be added to obtain the complete expression for WxAyAij ^xAy'Aki- 
The corresponding calculations are straightforward and the result reads 


( ( tiWx AyAy'A 

Jx±^ JA 


fid 


R 




Qdl 


,-o-Ry 


o-°-rv ' 
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(30) 

(31) 















where an = n'yCR/2 plays the role of a saturation scale. 

The above calculation, performed at large and for the case of a large target 
subdivided into many uncorrelated regions, has no immediate application to realistic ex¬ 
periments. However, it provides us with a set of non-perturbative inclusive and diffractive 
parton distributions which are highly constrained with respect to each other. Assuming 
that some of the essential features of diffractive and inclusive DIS are common to both 
the above model and the realistic proton case, we use the basic formulae Eq. (|30D and 


(1^ for a phenomenological analysis. For this purpose, H and a = n'yNc/4: are consid¬ 
ered as new fundamental parameters, giving rise to the following formulae for the basic 
hadronic quantities required in Sections ^ and |^, 


lx± J A 






1 - 






1 — e 


-2ay'- 


. (32) 

(33) 

(34) 


A similar, Glauber type y‘^ dependence has been recently used in the DIS analysis of . 
Note that according to Eqs. (|32|)-(p^ the diffractive structure function is not suppressed 
by a colour factor relative to the inclusive structure function, as originally suggested 
in 1^ 


5 Numerical analysis 

The model of the previous section provides an example for the relation between diffractive 
and inclusive parton distributions in the semiclassical approach. Although the derivation 
was based on a large hadronic target (with radius much greater than 1/A), we expect 
some qualitative features of the resulting distributions to apply to the proton as well. 
The above model distributions are used as non-perturbative input at some small scale 

and are evolved to higher using the leading-order DGLAP equations [^. The 
non-perturbative parameters of the model, as well as the scale Q^, are then determined 
from a combined analysis of experimental data on inclusive and diffractive structure 
functions. 

At hrst sight, the semiclassical description of parton distribution functions always 
predicts an energy dependence corresponding to a classical bremsstrahlung spectrum: 
q{x),g{x) 1/x. However, this naive prediction assumes the averaging over the soft 

field configurations inside the proton to be independent of the energy of the quark pair 
used to probe these configurations. As already outlined in Sect. ||, one expects that, in a 
more complete treatment, a non-trivial energy dependence is induced since the averaging 
procedure encompasses more and more modes of the proton field with increasing energy of 
the probe. We are, however, unable to calculate this non-perturbative energy dependence 
from first principles. Instead, we choose to parametrize it in the form of a soft, logarithmic 
growth of the normalization of diffractive and inclusive parton distributions with the 
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collision energy ~ 1/x, consistent with the unitarity bound. This introduces one further 
parameter, L, into the model: 

—>■ (L — Inx)^ . (35) 


Including this energy dependence, our model yields the following compact expressions 
for the semiclassical inclusive parton distributions Eqs. (|),(H) at a low scale Q^: 


xq{x,Qo) = 


xg{x,Ql) = 


2 , aQ,Nc {L — \n xY f Qq 


2aVtNc (L — Inx)^ 


In ^ - 0.6424 
a , 


(36) 

(37) 


Being derived in the semiclassical approach, these expressions are only valid in the small- 
X region, which we dehne by x < 0.01. In our numerical analysis, we shall multiply the 
above expressions with (1 —x) to ensure vanishing of the distributions in the limit x —>■ 1, 
which is required for the numerical stability of the DGLAP evolution. 

The corresponding expressions for the diffractive distributions Eqs. dl),® can be 
derived in a similar manner. The integrations over the momentum variables are outlined 
in Appendix A. One is then left with the expressions 


dq (/3,C,Qo) 
di 

dgiY,^,Ql) 


aUNYl-Y) (L-lnO%,,, 

27r3^2 hKP) > 

(38) 


(39) 


with the functions fq,g{Y) being dehned in Appendix A. The (3 spectrum of the diffrac¬ 
tive parton distributions at Qq is independent of the unknown non-perturbative param¬ 
eters. Note that our model does not specify whether, in the diffractive case, the energy- 
dependent logarithm should be a function of x or of Y However, both prescriptions differ 
only by terms proportional to In /3, which can be disregarded in comparison with In x or 
In^ in the small-x region. Like their inclusive counterparts, these distributions are only 
valid in the region ^ < 0.01, where we believe the semiclassical approach to apply. 

The above equations summarize our input distributions, depending on a, H, L, and 
on the scale Qq, at which these distributions will be used as boundary condition for 
the leading-order DGLAP evolution. At this order, the measured structure function F 2 
coincides with the transverse structure function discussed in Sects. |] and |^. In dehning 
structure functions and parton distributions, we assume all three light quark flavours to 
yield the same contribution, such that the singlet quark distribution is simply six times 
the quark distribution dehned above, both in the inclusive and in the diffractive case: 


S(x,Q2) = 6g(x,Q2), 


dE(e,/3,Q^) dg(e,AQ^) 

di di 


(40) 


Valence quark contributions are absent in the semiclassical approach, which does not 
account for the exchange of havour quantum numbers between the proton and the fast 
moving virtual photon state. 
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Charm quarks are treated entirely as massive quarks in the hxed flavour number 
which is appropriate since, in the range under consideration, charm 
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threshold effects are far more important than the resummation of InQ^/m^ terms. We 
thus fix nj = 3 in the DGLAP splitting functions and evolve only gluon and singlet 
quark distribution. The structure functions F 2 and are then given by the singlet 

quark distribution and a massive charm quark contribution due to boson-gluon fusion. 
Explicit formulae can, for example, be found in |^. For our numerical studies we use 
ALo,n /=3 = 144 MeV {as{Mz) = 0.118), rric = 1.5 GeV, irtb = 4.5 GeV, and we evaluate 
the massive charm quark contribution for a renormalization and factorization scale jjLc = 
2mc- 


The resulting structure functions can be compared with HERA data on the 
inclusive structure function F 2 {x,Q‘^) [Q and on the diffractive structure function 
[!]• These data sets from the HI and ZEUS experiments are used to 
determine the unknown parameters of our model. We apply the following selection crite¬ 
ria to the data: x < 0.01 and ^ < 0.01 are needed to justify the semiclassical description 
of the proton colour field; with Qq being a fit parameter, we demand a sufficiently large 

minimum = 2 GeV^ to avoid that the data selection is influenced by the current value 

of Ql; Anally we require > 4 GeV^ in the diffractive case to justify the leading-twist 
analysis. 

We determine the optimum set of model parameters from a minimization of the total 
(based on statistical errors only) of the selected data, using the MINUIT package 
The resulting set of parameters is 

Ql = 1.23 GeV^ , 

L = 8.16, 

G = (712 MeV)-2 , 
a = (74.5 MeV)^ . (41) 



The distributions obtained with these fitted parameters yield a good qualitative descrip¬ 
tion of all data on inclusive and diffractive DIS at small x, as illustrated in Figs. 1,1 
and All parameters are given with a precision which allows to reproduce the plots, but 
which is inappropriate with respect to the crudeness of the model. The starting scale 
is in the region where one would expect the transition between perturbative and non- 
perturbative dynamics to take place; the two other dimensionful parameters and a 
are both of the order of typical hadronic scales. 


Our approach fails to reproduce the data on F^^^^ for low (open dots in Figs. ^ 
and^). This might indicate the importance of higher twist contributions in this region, as 
suggested in p. It is interesting to note that a breakdown of the leading twist description 
is also observed for inclusive structure functions , where it occurs for similar invariant 
hadronic masses, namely < 4 GeV^. 


The perturbative evolution of inclusive and diffractive structure functions is driven 
by the gluon distribution, which is considerably larger than the singlet quark distribution 
in both cases. The ratio of the inclusive singlet quark and gluon distributions can be read 
off from Eqs. (P6[) and (^). With the fit parameters obtained above, it turns out that the 
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inclusive gluon distribution is about twice as large as the singlet quark distribution. In 
contrast, the relative magnitude and the (3 dependence of the diffractive distributions are 
completely independent of the model parameters. Moreover, their absolute normalization 
is, up to the slowly varying factor l/asiQl), closely tied to the normalization of the 
inclusive gluon distribution. 


Figure displays the diffractive distributions (multiplied by [3 and thus reflecting 
the distribution of momentum carried by the partons) for hxed ^ = 0.003 and different 
values of Q^. The [3 dependences of the quark and the gluon distribution at are 
substantially different: the quark distribution pdT,/d^ is peaked around f3 ~ 0.65, thus 
being harder than the distribution /?(! — /?) suggested in [^. It vanishes like /? for ^ 0 
and like (1—/?) at large (3] the gluon distribution pdg/d^, on the other hand, approaches a 
constant for [3^0 and falls off like (1 —/d)^ at large (3. This asymptotic behaviour in the 
small- and large-/3 region is in agreement with the results obtained in the perturbative 
approach of |^. In spite of the (1 — behaviour, gluons remain important even at 
large f3, simply due to the large total normalization of this distribution (the (3 integral 
over f3dg/d^ at Ql is approximately three times the {3 integral over j3dT,/d^). As a result, 
the quark distribution does not change with increasing for (3 ~ 0.5 and is only slowly 
decreasing for larger values of /?. 


The dependence of the diffractive structure function on (3 and is illustrated in 
Fig. H, where we compare our predictions with data from the HI and ZEUS experi¬ 
ments 0 at fixed ^ = 0.003 (HI) and ^ = 0.0042 (ZEUS). It must be pointed out that 
these data points are based on a combination of data taken at various values of ^ (in¬ 
cluding ^ > 0.01), which have been extrapolated to fixed The energy dependence of 
the diffractive structure function used for this extrapolation is different from the energy 
dependence employed in our model, so that a detailed comparison of the data with our 
results should only be made with some caution. Disregarding the large-/? region, our 
model gives a good description of the [3 dependence of the diffractive structure function 
for all values oi Q^. 


The validity of our approach can be tested by studying DIS observables other than 
the diffractive structure function used to fit our model parameters. In particular, it 
would be of interest to study quantities which are directly proportional to the diffractive 
gluon distribution, such as the charm content of The dashed lines in Fig. ^ show 

our prediction for the diffractive charm structure function due to photon-gluon fusion. 
Moreover, this structure function is predicted to have the same, non-perturbative ^ de¬ 
pendence as (3,Q‘^). The charm contribution is sizeable in the small-/? region. 

Since the full diffractive cross section at hxed x is obtained after integration with the 
measure d/?//?, this region yields a substantial diffractive charm cross section. 


The (3 dependence of diffractive parton distributions used in the present numerical 
analysis is a direct result of our specihc model for the colour held averaging. It would 
therefore be very interesting to repeat the analysis with more sophisticated models, such 


as the stochastic vacuum model that was utilized for dihractive meson production in [23 
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6 Conclusions 


In the target rest frame, DIS at small x can be viewed as the interaction of a partonic 
fluctuation of the virtual photon with the proton colour held. The semiclassical approach 
assumes this proton held to be dominated by soft modes, which have the sole ehect of 
introducing a non-Abelian phase factor for each parton. In this approach, a very similar 
description of inclusive and dihractive events emerges, the latter being realized if the 
partonic huctuation leaves the proton in a colour singlet state. Inclusive and dihractive 
structure functions are therefore calculated by evaluating diherent colour contractions of 
the same soft scattering amplitudes. 

Matching the semiclassical and partonic description of structure functions at some 
low scale Qq, where logarithmic corrections are still small, a set of semiclassical parton 
distribution functions for inclusive and dihractive processes is dehned. The evolution to 
higher scales is then determined by the conventional leading-order DGLAP equations. 
Initial quark and gluon distributions are expressed in terms of averages over the proton 
colour held. 

In the semiclassical framework, a very special role is played by the inclusive gluon 
distribution. In contrast to both the inclusive quark distribution and the dihractive quark 
and gluon distributions, it is only sensitive to the short distance structure of the proton 
held, and it is enhanced by an explicit factor l/ag. As a result, the observed dominance 
of the inclusive over the dihractive DIS cross section emerges. 

To study the semiclassical distributions in more detail, we introduce a non- 
perturbative model for the proton colour held. This model is derived for a very large 
hadronic target, which can be subdivided into diherent zones of uncorrelated colour held 
strengths. A non-trivial example of a non-perturbative set of semiclassical distributions 
results. In spite of its lacking theoretical justihcation in the case of a proton target, this 
set of distributions can serve as a basis of a phenomenological analysis. Our model de¬ 
pends on two free parameters which are related to the average held strength in a zone 
and to the total geometrical size of the target. A further parameter has to be introduced 
to account for the unknown energy dependence induced by the averaging over the held 
conhgurations. Thus, our model describes all semiclassical parton distributions by three 
unknown non-perturbative parameters and the matching scale Q^. 

These three model parameters and Qq are determined from a combined ht to mea¬ 
surements of inclusive and dihractive |Q] structure functions at small x. The resulting 
distributions yield a satisfactory description of the structure functions F 2 {x,Q^) and 
It turns out that the evolution of both structure functions is mainly 
driven by large gluon distributions. 

The good agreement of our model with the experimental data allows us to conclude 
that both inclusive and dihractive DIS at small x can be described in a unihed picture 
in the semiclassical approach. A simple model for the colour held averaging, adopted 
from a large target, is used to compute the semiclassical parton distributions at some 
initial scale Qq. The behaviour above Qq is then determined by perturbative evolution 
at leading order in and leading twist. 
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The rise of F 2 {x, Q‘^) and of /3, at small x have the same, non-perturbative 

origin in the energy dependence of the average over soft held conhgnrations in the pro¬ 
ton. With increasing Q^, this rise is enhanced by pertnrbative evolntion in the case of 
the inclnsive strnctnre fnnction, while it remains nnchanged in the diffractive strnctnre 
function. 

The obtained parton distribntions can be nsed to predict a broad spectrnm of ob¬ 
servables in diffractive and inclnsive processes nsing standard methods of pertnrbative 
QCD. The qnalitative relation between diffractive and inclnsive DIS and the nniversal 
^ dependence of large-mass diffraction are gennine predictions of the semiclassical ap¬ 
proach. By contrast, the /3 dependence of the diffractive distribntions is a resnlt of onr 
specihc model for the colonr held averaging. Crncial tests of this model can be performed 
in fntnre measnrements of dihractive hnal states. At the same time, other models for the 
averaging procednre can be tested in the present framework. 

We wonld like to thank M.F. McDermott, B.R. Webber and H. Weigert for valnable 
discnssions and comments. 


Appendix A: Diffractive parton distributions 


In the following, we collect several formnlae which are nsefnl in connection with the 
semiclassical expressions for the dihractive parton distribntions (Sect. ^ and with their 
explicit evalnation within onr model of the proton colonr held (Sect. 

In Ref. PI, the dihractive qnark distribution was given in the form 


dq{b,0 _ 2 rd^k'^k'^ r 

d^ ~ ^2(1-6)2 J (27r)8W -ix 


d?k± kx_ 

uk'"^ + 





2 




(42) 


where u = h/{l — h). Starting from this result, Eq. (|T^ of Sect. || can be derived using 
the identity 


r d?k± 

J (27r)2 m + kl 


2vr y 


NK,{yN) . 


(43) 


For the phenomenological analysis, we use the specihc model of Sect. Inserting for 
tTW^^{y±^)tTW^}{y'^) in Eq. (|42|) the expression from Eq. (^), the y± integration and 
some of the momentum integrations can be carried out. The result reads 


dqjb,^) _ anNcjl - b) 
di 27r3e2 ^ ’ 


where fq{h) is an integral over two Feynman-type parameters, 

( _ Vh + X _\ ( _ Vh + x' 


f,(h) = 4 hi-t+(v^+x)y 


(44) 


(45) 
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We were not able to obtain an analytical expression for this integral. It is, however, easily 
evaluated at 6 = 0 and 6 = 1 yielding fq{0) = 1/2 and /g(l) = Stt^/S — 2 . 

Similar formulae hold for the diffractive gluon distribution. The result of reads 


dg{b,0 ^ b r r 


d?k± tij 

uk'‘^ + k‘^ 


'y± 




(46) 


where 


Using 


tij dij 2 


kik 




uk'"^ 


f d^k±_ (^ ^ ^kikj\ 

J (27r)2 + ) N^ + kl 




(47) 

(48) 


one obtains Eq. (0) in Sect. ^ Inserting for trlU^(?/_L)trhU^l(?/^) the expression in 
Eq. (^41) of Sect. the y± integration and some of the momentum integrations can be 
carried out. The result has the same structure as the diffractive quark distribution. 


dg{b,0 

d^ 271^^% ’ 


(49) 


where fg{b) is given by the two-dimensional integral 


m 



dxdx' 


( 1 — 6-1- 3(1 -I- x^b \ / 1 — 6-1- 3(1 -f- x^b 

y (1 -t- xY{l - 6 (1 xYbY ) \{l + - b+ {1 + x'Ybf 

{x + xjh + (1 - 6 ) (y^ + Y+l?) 


(50) 

This integral is easily evaluated for 6 = 0 and 6=1 yielding /g(0) = 41n2 and fg{l) = 
457r^/32 — 17/2. For general 6, we have evaluated /g(6) and fg{b) only numerically, the 
results can be inferred from the solid curves in Fig. 


Appendix B: Comparison with a perturbative model 


It is the purpose of this appendix to outline how the perturbative results of can 
be derived on the basis of the semiclassical formulae for diffractive parton distributions 
of m (cf. Eqs. (^ and (||) of this paper). 

The authors of [Q study diffraction as quasi-elastic scattering off a special target 
photon that couples to only one flavour of very massive (M 3> A) quarks. The large 
quark mass justihes a completely perturbative treatment of the target and the diffractive 
system. In this situation, the required t channel colour singlet exchange is realized by two 
gluons coupling to the massive quark loop of the target. In the semiclassical approach, 
these two gluons are understood to be radiated by the massive quark loop and are treated 
as the colour held generating trfUtrlUh The semiclassical calculation proceeds as follows. 
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(51) 


Equations (|4^ and (|^) have the structure 


dfa 



trhh,^ trhhj^ 


where Fa (with a = q, g) is a. linear functional depending on /trlEEj^(j/^) trlEl^(?/(i_), 
interpreted as a function of y± and y'j_. To be differential in t, one simply writes 


dfa 

dt 







(52) 


with = —t. 

The held responsible for trhT is created by a small colour dipole which, in turn, is 
created by the special photon that models the target. At leading order in perturbation 
theory, the colour held of a static quark is analogous to its electrostatic Coulomb held. 
The held of a quark travelling on the light cone in X- direction (x± = xq ± ^ 3 ) at 
transverse position 0 _l has therefore the following line integral along the x+ direction. 


2 


A_ dx+ = —ig^ 


d'^k± g*^-L^-L 

(27r)2 kf 


(53) 


It is exactly this type of line integral that appears in the exponents of the non-Abelian 
phase factors U and W that form W (cf. Eq. (3) of |^). A straightforward calculation 
shows that the function trW produced by a dipole consisting of a quark at p± and an 
antiquark at 0 j_ reads 


tTW^^{y±) 


g\N^ - 1 )Tr \ f (1 - (l - 

2 {2TTYk\k'‘l 


(54) 


where Tp = 1/2 and Ta = have to be used for the fundamental and adjoint represen¬ 
tation respectively. 

The hnal formulae for the dihractive parton distributions of the target are obtained 
after integrating over the transverse sizes of the colour dipoles with a weight given by 
the qq wave functions of the incoming and outgoing target photon. They read 


dfa 

d^ dt 


= J dz d p± J dz d p^ 


'X± Jx , 


trhl4^ trirj, 


(55) 


^ i H d±) 'f-f{z,p±,0x,e±)] pf,pf, ef) Op, ep) 


where tiWx^iyx) is produced by the held of a quark at p± and an antiquark at Op, and 
trWx' {y'jJ is produced by the held of a quark at pf and an antiquark at Op, as detailed 
in Eq. (|5j[) . 
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The wave function '^ 7 ( 2 ;, pn, 0 _l, ej,) characterizes the amplitude for the fluctuation 
of the incoming target photon with polarization e and transverse momentum 0 ^ into a 
qq pair with momentum fractions 2 ; and 1 — z and transverse separation p^. Similarly, 
the wave function characterizes the amplitude for the recombination of 

this qq pair into a photon with polarization e' and transverse momentum = —q±. The 
summation over the helicities of the intermediate quark states, which are conserved by 
the high-energy gluonic interaction, is implicit. 

The required product of photon wave functions can be calculated following the lines 
of 0,|4|. It reads explicitly 


e'^) ^^7(2;, p±. Ox, ex) 


(56) 


Nce‘^el 

2(27r)5 


/ tr<ht(x, fcl, M, e'^)$(x, fcx, M, ^ 

Jk I ,k' 


where we have used the notation of ilH 


^{z,k±,M, ex) 


1 

{kl + M2) 


[{1 — z) e± ■ a k_}^ ■ a — zkj^ ■ a ■ a + iM €± ■ a] , 


(57) 


M is the quark mass and ai ^2 are the hrst two Pauli matrices. Note that for p'j_ = 0, 
the average of the diagonal elements (ex = e'j^) in Eq. (|56D reproduces the well known 
formula for the square of the photon wave function Q. 

Inserting Eq. (^) into Eq. ( ^5]) and introducing explicitly the required functionals 

0 


Fa specihed by Eqs. 


and 


the formulae of 


for diffractive quark and gluon 


distribution are exactly reproduced. 


References 


[ 1 ] HI Collab., S. Aid et al., Nucl. Phys. B470 (1996) 3; 

ZEUS Collab., M. Derrick et ah, Z. Phys. C72 (1996) 399 

[2] HI Collab., C. Adloff et ah, Z. Phys. C76 (1997) 613; 

ZEUS Collab., J. Breitweg et ah, preprint DESY-98-084 (Piep-ex/9807010|) 

[3] N.N. Nikolaev and B.G. Zakharov, Z. Phys. C49 (1991) 607 

[4] K. Golec-Biernat and J. Kwiecihski, Phys. Lett. B353 (1995) 329; 

A. Edin, G. Ingelman and J. Rathsman, Z. Phys. G75 (1997) 57; 

E. Gotsman, E. Levin and U. Maor, Nucl. Phys. B493 (1997) 354; 

A. Bialas, R. Peschanski and G. Royon, Phys. Rev. D57 (1998) 6899 


[5] 

[ 6 ] 


T. Gehrmann and W.J. Stirling, Z. Phys. G70 (1996) 89 

J. Bartels, J. Ellis, H. Kowalski and M. Wiisthoff, preprint GERN-TH-98-67, DESY 
98-034 and DTP-98-02 (|hep-ph/9803T97D 


20 












[7] W. Buchmiiller, M.F. McDermott and A. Hebecker, Nucl. Phys. B487 (1997) 283; 
ihid. B500 (1997) 621 (E) 

[8] A. Hebecker, Nucl. Phys. B505 (1997) 349 

[9] G. Ingelman and P.E. Schlein, Phys. Lett. B152 (1985) 256; 

L. Trentadue and G. Veneziano, Phys. Lett. B323 (1994) 201; 

A. Berera and D.E. Soper, Phys. Rev D50 (1994) 4328 

[10] V.N. Gribov and L.N. Lipatov, Sov. J. Nucl. Phys. 15 (1972) 438, 675; 

G. Altarelli and G. Parisi, Nucl. Phys. B126 (1977) 298; 

Yu.L. Dokshitzer, Sov. Phys. JETP 46 (1977) 641 

[11] L. McLerran and R. Venugopalan, Phys. Rev. D49 (1994) 2233 

[12] A. Hebecker and H. Weigert, Phys. Lett. B432 (1998) 215 

[13] W. Buchmiiller, Phys. Lett. B353 (1995) 335; 

W. Buchmuller and A. Hebecker, Nucl. Phys. B476 (1996) 203 

[14] F. Hautmann, Z. Kunszt and D.E. Soper, preprint GERN-TH-98-168, ETH-TH/98- 
14 and GITS 648 (|hep-ph/9806298|) 

[15] J.G. Gollins, Phys. Rev. D57 (1998) 3051 

[16] J. Jalilian-Marian, A. Kovner, L. McLerran and H. Weigert, Phys. Rev. D55 (1997) 
5414 

[17] K. Golec-Biernat and M. Wiisthoff, preprint DTP/98/50 (|hep-ph/9807513| ) 

[18] W. Buc hm iiller and A. Hebecker, Phys. Lett. B355 (1995) 573 

[19] M. Gliick, E. Reya and M. Stratmann, Nucl. Phys. B422 (1994) 37 

[20] F. James, MINUIT, GERN Program Library Long Writeup D506 (1994) 

[21] A. Milsztajn and M. Virchaux, Phys. Lett. B274 (1992) 221; 

A.D. Martin, R.G. Roberts, W.J. Stirling and R.S. Thorne, preprint DTP-98-52 
(|hep-ph/9808371|) 

[22] A. Donnachie and P.V. Landshoff, Phys. Lett. B191 (1987) 309; ibid. B198 (1987) 
590 (E) 

[23] H.G. Dosch, O. Nachtmann and M. Rueter, preprint HD-THEP-98-22 and TAUP- 
2502-98 (|hep-ph/98063421) 

[24] M. Wiisthoff, Phys. Rev. D56 (1997) 4311 


21 










1.5 

1 

0.5 

0 

1.5 

1 

0.5 

0 

1.5 

1 

0.5 

0 

1.5 

1 

0.5 

0 

1.5 

1 

0.5 

0 


Figure The inclusive structure function F 2 {x,Q‘^) at small x computed in the semi- 
classical approach, using the fitted parameters given in the text. Data taken from |^. 
The data with = 1.5 GeV^ are not included in the £t. 
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Figure The diffractive structure function at small ^ computed in the 

semiclassical approach, using the htted parameters given in the text. HI data taken 
from 0], The open data points correspond to < 4 GeV^ and are not included in the 
£t. 
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Figure The diffractive structure function at small ^ computed in the 

semiclassical approach, using the fitted parameters given in the text. ZEUS data taken 
from 0. The open data points correspond to < 4 GeV^ and are not included in the 
£t. 
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Figure ^ Diffractive quark and gluon distributions at the initial scale and after 
evolution. 
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Figure |^: Dependence of the diffractive structnre fnnction on {3 and compared 
to data from HI (left) and ZEUS (right) p. Open data points correspond to < 

4 GeV^. The charm content of the structure function is indicated as a dashed line. 
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